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We consider a generic time-reversal invariant model of fermions hopping randomly on a 
square lattice. By means of the conventional replica-trick within the fermionic path-integral 
formalism, the model is mapped onto a non-linear a-model with fields spanning the coset 
U(4N) I Sp(2N), N — > 0. We determine the proper scaling combinations of an infinite family 
of relevant operators which control deviations from perfect two-sublattice symmetry. This 
allows us to extract the low-energy behavior of the density of states, which agrees with earlier 
results obtained in particular two-sublattice models with Dirac-like single-particle dispersion. 
The agreement proves the efficacy of the conventional fermionic-path-integral approach to 
disordered systems, which, in spite of many controversial aspects, like the zero- replica limit, 
remains one of the more versatile theoretical tool to deal with disordered electrons. 
PACS numbers: 73.20. Jc; 73.20.Fz; 71.30.+h 

I. INTRODUCTION 



It is known that localization does not occur in any dimension at the band-center en ergy 
of tight binding models on bipartite lattice-Hamiltonians 

QflflflBBBBBH ill. 

Gade and Wegner 4 , 0] first realized that these models correspond to a particular class of 
non-linear a- models in the zero replica limit, so called two-sublattice models, which exhibit 
an additions! chiral Symm e try fl £ Q . As a eonse q nence of this s ymm e try , they were able 
to prove that, when the chemical potential is right in the centre of the band, quantum- 
interference corrections to the /3-function vanish exactly. In addition they showed that, 
unlike conventional disordered systems, the density of states (DOS) p(E) near the band- 
center E = is strongly affected by disorder. They actually predicted a diverging behavior 



Q, (see also Refs. Q, Q ) 

where A and B are positive constants. 

The subleading dependence exp ( — Ayjhx \B/E\) has been recently questioned by 
Motrunich, Damle and Huse Q|. They analysed the strong disorder limit of a set of models 
which belongs to the two-sublattice class and found that the correct subleading dependence 
is instead exp ( — A\ InB /E\ 2 ^ 3 ). This result was later confirmed by field-theoretical ap- 

n n 

proaches based on supersymmetry |13|] and on replica trick applied to the so-called 
Hatsugai-Wen-Kohmoto (HWK) model, which describes electrons hopping randomly on a 
square lattice in the presence of a 7r-flux per plaquette. This model is particularly suitable 
for a weak-disorder field-theoretical approach. Indeed, for uniform hopping, the low-energy 
single-particle spectrum of the HWK model is composed by two Dirac-like cones, which al- 
lows the use of the full machinery of Conformal Field Theory when a weak random-hopping 
component is included. The important breakthrough put forward by these analyses is that 
many disorder- average quantities, like the density of states, are determined by an infinite 
set of relevant local operators with negative dimensions 13] . 

In reality, these new results raise an intriguing question about the concept of universal- 
ity commonly accepted in disordered systems, according to which the HWK model should 
be representative of any two-sublattice model since the action of its long-wavelength diffu- 
sive modes is a non-linear a-model in the same universality class as g eneric two-sublattice 
models. However the results found by Mudry, Ryu, and Furusaki [l3j] and by Yamada and 
Fukui have been obtained working directly with the HWK Hamiltonian without going 
through the non-linear a-model mapping, just thanks to the Dirac-like dispersion. Lacking 
an independent derivation starting from the non-linear cr-model, it is not a priori obvious to 
what extent these results are actually generic to any two-sublattice model. 

In this work, we are going to show that the results of Refs. |3, [3, [l^ can be fully 
recovered through the conventional non-linear cx-model approach based on the replica-trick 
within the fermionic path-integral formalism, without assuming any Dirac-like dispersion. 
Besides satisfying a purely theoretical curiosity, this result proves once again the strength 
of the conventional approach to disordered systems. 

The paper is organized as follows. In Section II we introduce the model as well as the non- 
linear a-model which describes its long-wavelength diffuse modes. In Section III we analyse 
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the scaling behavior of some operators which are compositions of such diffuse modes, called 
soft operators, and whose /3-functions are solved in Section IV. Finally in Section V we apply 
what found in the previous sections to calculate the mean density of states near the center 
of the band. 



II. THE MODEL 



A two-sublattice model is described by a non-interacting Hamiltonian on a bipartite 
lattice, with sublattices A and B, of the general form 

H = -J2 J2 tRR,c R c R' + H - c - C 1 ) 

ReA R'eB 

with random matrix elements tRR> which only connect one sublattice to the other. It follows 
that, if ^(R) is an eigenstate with eigenvalue E, then the wavefunction <p{R) ^{R), where 
4>(R) = 1 if R G A and 4>{R) = — 1 if R G B, is also an eigenstate with eigenvalue —E. This 
also implies that any eigenstate at E = is doubly degenerate, unless boundary conditions 
break the degeneracy between the two sublattices. Following Ref. Q], we introduce, within 
the fermionic path-integral formalism, the following Nambu spinors 




where cr and cr are Grassmann variables with components CR tPta and cr iPi(1 in which R refers 
to a lattice site, p = ± is the index of positive or negative frequency components, a — 1, .., N 
is the replica index which has to be sent to zero at the end of the calculation. In addition 
we define conjugate fields through 

^R = [c^r] 1 = - Cr,Cr), 

where c is the charge conjugation matrix c = —1T2. Here and what follows the Tj's, 
i = 1, . . . , 3, are Pauli matrices acting in the Nambu space. With these fields the ac- 
tion describing the model at a fixed disorder configuration for fields with energy E ± iuo is 

S = - ^2 i>R{ ES RR' ~ ius 3 5 R R, - tRR')ipR, (2) 
ReA R'eB 
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where S3 is the third Pauli matrix acting on the two frequency components E ±iu and the 
frequency transferred 2a; plays the role of a symmetry breaking field. Since the detailed 
derivation of the non-linear a-model starting from the action Q is known, see for instance 
Ref. 0, here we just outline the main steps emphasizing the peculiarity of two-sublattice 
models. We start noticing that the action when E = uj = is invariant under the transfor- 
mation 

which is allowed within the path-integral formalism since the Grassmann variables cr and cr 
are independent. This additional abelian symmetry plays an important role in these models, 
as originally recognized by Gade and Wegner |4f. Within the replica-trick technique, the 
average of disorder can be performed directly on the action (j2J) and generates a non-local 
interaction which connects two sites belonging to different sublattices. This interaction is 
then decoupled by an Hubbard- Stratonovich transformation, introducing an auxiliary field 
Q. Due to the non-locality of the interaction, in the long-wavelength limit this field has 
two components, one uniform, Qo(R), and the other staggered, 4>{R)Q^{R), where both 
Qo and Q3 vary smoothly in space. In particular, if we introduce the Pauli matrices ji, 
i = 1,2,3, as well as the identity matrix 70 in the two-sublattice space, the action of 
73 being just that of <f>(R), the auxiliary field to which the electron density is coupled is 
Q(R) = Qo(R) 7o + iQz{R) 73- Both Q and Q 3 are AN x AN hermitian matrices in Nambu, 
energy and replica spaces [AN = 2 (Nambu spinor dimension) x 2 (frequency components) 
xN (replicas)]. The derivation proceeds then through the following steps: i) integrating over 
the Grassmann variables, ii) expanding the effective action around the symmetry breaking 
saddle point Q = £33, where S is the inverse relaxation time, iii) integrating out massive 
modes and just focusing on low-energy long-wavelength transverse fluctuations. In this way 
one obtains the following effective action for the auxiliary field when E = [5], 0, 



F = I dr 



^Tr (VQt( r )VQ(r)) - ( Tr {QHr)VQ(r) l3 }) 2 - ^ Tr (s 3 Q(r)) 



(3) 

In the above equation we have rescaled Q — > SQ with Q = W S3U and the unitary transfor- 
mation U E U(AN) /Sp(2N) fulfilling 

^ t = 7 1 C/ t 7 1 = cC/*c*. (4) 
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The latter relation implies the constrain Q^Q = 1. The coupling a is the conductivity in 
Born approximation, II is related to the staggered density of states fluctuations and v 
is the density of states at the chemical potential. The transformation U can be written as 
U = e w with W = W / °7o + W /3 73. The charge conjugation invariance implies for the Q-fields 
the following relation 

cQ l & = Q. (5) 



A. Gaussian propagators 

In the Wilson-Polyakov renormalization group (RG) approach one assumes a separa- 
ble form for the transformation U = U(j\U( 8 \ where U(j\ involves fast modes with momentum 
q E [A/s, A], while Ur s > involves slow modes with momentum q G [0, A/s], with A being the 
high momentum cut-off, and s > 1 the rescaling factor. By introducing upper indices, 
a,b,c, in Nambu space which assume the values 1 and 2, and lower multindices in replica 
and energy spaces (n = (r,p), with r being the replica index and p the sign of the frequency) 
we expand Q to the second order in VK-fast obtaining 

r^ad ^ r^ad _i_ Tj\ab q T/T/&C TT C< ^ _l_ ^_TT^ a ^ Q T/f/^ c ' T/f/ c ' c TT cf ^ (P,\ 

Vnin 2 — ^riYri2 (s) ' U nim 1 "mi VV m 1 m 2 U m 2 n 2 ' n H">1 mi Hi! N Hj^ ' V / 

In the above equation the sums over repeated indices are assumed, the label (s) is dropped 
in U and S n — p is the sign of the infinitely small frequency related to the index n. The 
gaussian propagator reads 

(W:£ b (k)W« p cd (k)) = ((-) a S n S m - 1) D a (k) {5 al 5 bl (5 cl 5 dl 5 np 5 ml 

(SnpSmi + (-TSniS mp ) + 5 a2 5 b2 (5 c2 5 d2 5 np S ml - {-) a S cl 5 dl 8 nl 6 mp ) 

+T5 a3 5 ab 5 cd 5 nm 5i p /2} , (7) 

where a = 0, 3 and 

D " {k) = D3{k) = ^F' (8) 

r =7TM' < 9 > 

with N are the number of replicas that we send to zero. The relation in Eq. (J3J) implies 

u£n = (* ab - 1) uit + 5 ab u^ lb+l = H a+b u^ la+ \ (io) 
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with the clock rules 1 + 1 = 2 and 2 + 1 = 1 for the upper indices. From Eq. (J7J) and Eq. 
(JTUJ) we obtain 

(Qt in2 ) = Qtn 2 + L^(2-8N + T) Q% n2 , (11) 

+L 3 [(2 -8N + 2Y)Qt n2 Qt ni + R b+C Q%£<j£H ~ Q%J&n*] > ( 12 ) 

where L a = J A , j0^D a (k) and the label (s) is dropped in all the Q-fields on the right-hand 
side of the equations. Eq. (JHJ) yields the properties 

Qnm Qmm (13) 

Qnm ~ ~Qmni (^) 

which imply some further properties, for instance, 

E( _\a+cf\ac+lf\a+lc \ * r\ac f\ca 

V J Vnin 3 Vn 2 n 4 / j ^ n\n?C^ nmii \ LO J 

ac ac 

useful to evaluate product of traces. From Eqs. (fT2j) and (JT5J), we get 

(J2QZn 2 Qn 3 n 4 ) = £ { ^ + 2L Q + HI <C*«!U. " ^ [Q^Qn^ + Q^Q^J } ,(16) 
ab ab 

(J2<n 4 Qnsn 2 ) = £ {[1 + 2L(1 + T)] Q^Qn^ + ^ [Q&jft^ ~ QZn 2 Qn 3 n 4 } } ,(17) 
ab ab 

(/] ( ?mn 3 ( ?n°n 2 ) = ^ {[1 + 2 ^(1 + H] Qn^mQn^ + 2 ^ [ ( 5mn 4 ( 5nln 2 _ Qm^M^nUJ } '(1 § ) 
ab ab 

where N = and L = Lq = L3. From the equations above we notice that the one- 
loop calculation leads to transpositions of the indices under the sum over Nambu space 
and that Eqs. (|16H18|) form a closed set of equations. Defining v\ = Y^abQnimQn3TH> 
V2 = ~ J2ab Qn^Qnani and v 3 = ~ Eafe Qtr^Qn^ the equations above can be summarized 
by 

(vi) = v i + 2LTv i + 2L(v 1 + v 2 + v 3 ). (19) 

Applying a rotation we can obtain three independent scaling operators V\ = v 1 + v 2 + V3, 
v 2 = V\ — v 2 and t>3 = v 2 — v$, with the following scaling behaviors 

(v 1 )=v 1 + 2L(3 + T)v 1 , (20) 
(v^ = Vi + 2LTvi , for i = 2, 3. (21) 
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From Eqs. and (|12[) one can evaluate the scaling behavior at one loop level of all the 
composite operators for the model © with symmetry U(AN) / Sp(2N). Another result which 
can be easily obtained from Eqs. (fTTj) and (|T2|) is the mean value of the full product of n 
matrices Q in the limit N — > 



(Q n > = Q n + L 



n-1 



— (2 + T)Q n -nJ2Q n ~ l Tr'Q l 



1=1 



(22) 



where Tr' is a trace which does not act on the sublattice space and T — H/a. 

In the absence of the sublattice (or chiral) symmetry, the transverse modes take values 
in the manifold Sp(2N)/ Sp(N) x Sp(N). In such a case the gaussian propagator ((7j) is valid 
only for a = 0, implying that (Q^n 2 ) = Qn\n 2 and 



ad ^icb 1 

n 3 n 2 \ ' 



while Eq. (122)) is replaced by 



(Q n )=Q n + L[-(n-l) 



n 



(23) 



(24) 



III. SOFT OPERATORS 



Let us consider now the following linear combination of moments of Q 



P « = E e c.nw 



i=l 



(25) 



where VW G {l,...,n}, all the sets of M positive ordered integers {rii} M = {ni, n^-} 
are such that Y2~iLi n i = n i namely {rii] M is a partition of n in J\f terms. The coefficients 
A <AA) are symmetric with respect to any transpositions of the indices, X (M) = 

({ni) M ) J r J r (n 1 ,..,nj,..,n i ,..,n Jt f) 

A (A ° . The number of couplings is given by the number of partitions of n, 

(n 1 ,..,n i ,..,nj,..,n A/ -) ^ 

2~^at=i S < n i>AA 1 = Pn- The trace Tr' is over all the degrees of freedom except those of 
the sublattice space in order not to miss operators induced by the renormalization. In this 
way, as we will see below, we have in Eq. ()25j) a complete set of operators transforming one 
to the other under the action of the renormalization group. 
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From Eqs. and ()12|) we get the following scaling behavior 



Af 
i=l 



1 + + 



N 



n Tr '^ r 



(26) 



Af 



1=1 
Af 



Af 



AJ2 n i n J Tr'Q n * +n > f[ ^'Q nk ) ^n Tr '^ nj J2 ^'Q^Tr'Q* 



I i>3 



i=l 



p=l 



for a generic product of traces of any power of Q. In the latter equation £) i>3 - * s a double 
sum and represents X^iX)! - =i- m the following we give two particular examples: i) For 
Af — 1, ni = n, using Eq. (|2b|). we get 



(Tr'Q") 



1 + L n + n 2 - 



n-l 



Tr'Q n - L n Tr'Q"- p Tr'Q ? 
p=i 



ii) For Af = n, then Vz n,; = 1, and we obtain in this case 

,r 



((Tr'Qr) 



1 + L I n + n 



(Tr'Q) n - 2Ln(n - 1) Tr Q 2 (Tr Q) n ' 2 



In the first term of Eq. (|26j) the product of the Af traces is reproduced, in the second 
term the number of traces instead is decreased by one, it is Af — 1, while in the third term 
the number of traces is increased by one, Af + 1. In any case the sum of the exponents of 
Q is equal to n. We can say, therefore, that for each positive integer n, the set of all the 
equations (|2*S|) . with all positive integers Af and {n^} such that ^2 i=1 rii = n, is a closed set 
of RG equations whose number is equal to the number p n of partitions of n. 

The corresponding p n /5-functions for the couplings, d^\ } ^/dlns, are the following 



(3 (ao 



A \ (A ° 

xvvj — a A, -+- q 



(f>' + »'§)C 



"A/0 



-4 K (M+1) A (AA+1) n(n, - n) 

i=l n=l 
A/" i-1 ^ 

— (1 — 5a/i) > > — n^ - \ ' (n» + n,-) 

\ ] ' / j / j J^(AT) (n 1 ,..rf.-...n.+n.-...n.rl V 1 J ' 



, 1 ^ 

= 1 J = l 



(27) 



where ni > n2 > ... > n^r and g = -^r^ is the resistivity. In Eq. (J2?]) the factor K n n is 
defined by 



K 



(Af) , (AQ A 
(AT) m ni l m n, "njrij/ 



2m 



(AT-1) 



(28) 
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where is the number of times the value of the integer rii is repeated in the string 

(ri\, ..,njsf) or equivalently, = (Im~ lm) where Im > l m £ N such that i G + / m } 
and V/ G {/a/ + 1, lm — 1} the /-cycle is absent in the class of permutations corresponding 
to {jii\ M , while m^ + ^ in the denominator is the multiplicity of rii + rij (the sum of the lower 
indices of K W ) in the string (m, ., fij, ., rij + rij, .,njv") that is equal to its multiplicity in 
the string (n 1; ..,n_\f) plus one. The same definition is valid for K_ n __, namely 

(./V+l) , (AT+1) . v 
(Af+i) _ ^-n \ m rn-n ~ n,ni-n) , . 

2m% } ' 1 J 

where r?7,^ +1> and rri^^ are the the multiplicities of the values n and — n respectively 
in the string (m, ..,n, — ft, .,nj^,rX) and m^. 5 in the denominator is the multiplicity of the 
value rii (the sum of the lower indices of K^]_ } _) in the string (ni, .., f4j /- fl, ., rij, ., rajy) that 
is the original string (n l5 ..,njv). 

Eq. f)26j) is equivalent to Eq. ()27|) since the first describes the scaling behavior of the operators 
while the second of the couplings. To show how Eq. (|27j) can be derived from Eq. let 
us consider the equations 

AT AT 

(Tr'Q n<+ ^ J] Tr'Q"*) = [...] - L (m + Uj ) \[Ti'Q n \ 
(Tr'Q n *- n Tr'Q n JJ Tr'Q^') = [...]- ALn(m - ft) JjTr'Q n \ 

*=1 

obtained directly applying Eq. (J2BJ). From these equations it is easy to see how the two cou- 
plings A^. i !n . + „. r)njv) and A^ ) nj) , >n .^ Mn) , corresponding to the two operators above, 
appear in the /3-function for A^' s, the coupling of Yl^Ti'Q^. 

The factors K^,^. take care of the symmetry of the couplings with respect to any trans- 
positions of the indices. The numerator in Eq. (}2*%j) is given by the number of different pairs 
(rii, rij) one can couple starting from objects of type rii and objects of type rij. If 
rii 7^ n ji the number of pairs is rn^rri^ while, if rij = rij, the number of pairs is given by 

(AO W i 

( m 2 ) = ^WT^ = l m l?( m l? ~ X )- The denominator in Eq. (J2Hj) is for rij = 

while 2m^. + ^. for rij 7^ n.,-, the factor 2 comes from exchanging rij <-» rij. 

By this procedure we can write down the one loop RG equations of the couplings of a 
generic product of traces of powers of the field Q. Now if we were interested to find one loop 
scaling operators it would be enough to find the real solutions of the p n — 1 independent 
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equations, algebraic through Eq. (|2*7j). among the equations 



A , P,(Af) — X. , (30) 

with the constrain Ylf n i = Ylf n i = n or alternatively to diagonalize the matrix M 
constructed by the coefficients of A's in the right-hand side of the equations (|27j) that has 
clearly rank p n . Indeed, denoting with A the p n -vector formed by the couplings A^j, Eq. 
(|27j) can be written in the following way 

x = (d + gM)X. (31) 

Calling A (A ° the real solutions of Eqs. (1301) (or alternatively the columns of the invertible 
matrix T that diagonalizes M), the resulting p n operators 

n Af 
A/=l {nj^ 1=1 

with z n = 1, ..,p n , are one loop scaling operators with the following scaling behaviors 

(O in ) - 



1+ L + ^n 2 +2j2n i (n i -i)j L O in , (33) 

where n\ > ri2 > ... ~> n^f is a partition of n. The factor in front of L is the z n th element of 
the diagonal matrix T _1 MT and goes from the value [(£ + 2) n 2 — n\, the most relevant, 
related to the partition {n}, to the value [(^ — l) + 2n] , related to the partition {1, 1}. 
Eq. (|3*3*|l is exactly what we can find also using the Young tableaux, already adopted, for 
instance, to evaluate the average of moments of the eigenf unctions of a particle in a random 
potential near the mobility edge [l?! ]. 

As in that case, one-loop renormalization leads to transpositions of indices of the matrix 
field under the sum of upper indices, taking advantage of the charge conjugation condition 
(0) (see, for instance, Eqs. (|16II18|) ). Indeed the factors n + 2 Yl'f n i( n i ~ *) are the 
eigenvalues of the operator X/»>j(v) w hich is the sum of all the transpositions (ij) acting 
on a set of n pairs {(1, 2), (3, 4), .., (2n — l,2n)}, the indices of n Q-matrices, and where 
(ij) is an identity on each pair. The dimension of the space spanned by this operator is 
(2n — 1)!! which is equal to the sum of the dimensions of some irreducible representations of 
the symmetric group S^n, those related only to the even partitions of 2n jl?) ]. The number of 
even partitions of 2n is also given by p n , the number of all the partitions of n. For instance, 
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if n = 2 the dimension of the space is 3, (the tree vectors Vi in Eq. (HHJ)) but the eigenvalues 
are two (see Eqs. (J2U1 12T|) . v 2 and t> 3 are degenerate). Indeed the irreducible representation 
related to the partition {4} has one dimension while the one related to the partition {2, 2} 
is a two-dimensional representation ["Lt1| . The additional term, \n 2 , results from the sum 
of two contributions, the first coming from the mean values ~(WWWU) , in all the terms 
(Q), which give \n and the second coming from (WWU...WWU), appearing in in all the 
average values (QQ), which give ^n(n — 1). 

Finally we can rewrite P n of Eq. ()25|) in terms of such scaling operators 

Pn 

P n = Y,^A. (34) 
For the sake of clarity we refer to an explicit example in Appendix B. 



IV. RG SOLUTIONS 

As we have seen before, we have decoupled all the moments of Q writing them in terms 
of the scaling operators Oi n , with i n = 1, ..,p n , which are coupled only to the equations of 
g, the resistivity, and T. In terms of the couplings ai n of the operators Oi n , the whole set of 
one-loop RG equations in e = d — 2 expansion Wi is the following 



Pg = -£9, 
Pr = 

Pa, n =da in + (A in + B n T)gai 



(35) 
(36) 
(37) 



where the coefficients Ai n and B n are defined in the following way 



Aj n = n + 2y j n i {n i - i), with ni > ... > n^, 



i=i 



B n = n 2 /2. 



(38) 
(39) 



These decoupled RG equations can be solved easily obtaining the following solution 



9 = Qos , 

r = r + —(i - s 



In 



a i n 0. 



d In s + — 

e 2 



9os 



-2e /e 



- 1) A in es e + B n (eT s e + 2g (s e - 1)) 



(40) 
(41) 

(42) 
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In 2-dimensions the solution of the RG equations is the limit e — > of Eqs. (j40H42[) . 



9 = 9o, 

r = T + 4^o In s, 



In 



a>i n o 



(43) 
(44) 
(45) 

(46) 



(2 + Ai n g + B n g T ) In s + 2B n g 2 (\n s) 2 . 

The flow equations for the original couplings in Eq. (|25p. called here Aj n with i r 
for simplicity, are 

Xi (s) = V T; j a, = s B " lns s A ^T~) A/ , 

»n V / / j l njn Jn / j i-njn Jn^n n ' 

where T is the invertible matrix that diagonalizes M and 

A tn =2 + A ln g + B n g T , (47) 
B n = 2B n g 2 . (48) 

Considering, in the limit s — > oo, only the first most relevant exponent, denoted by A Tn , for 
which a rn o 7^ 0, we have simply 

\ in (s) ~ S A - +BttIns T inr „ ^T-j n A, n (49) 
and if Aj n = A in 5 iniin we can write 



In 



KM 

A,- 



In 



.°r„0. 



(50) 



Inverting Eq. (|49|) we find that the couplings Aj n (s) reach an upper value A at the length 
scale 

1 



s = exp 



2 B r , 



, A r - AB n In I — — 



(51) 



Notice that in general terms the flow to strong coupling regime can be tuned and 
slowed choosing some particular starting configurations of the couplings, namely when 
V\ T~) Xi = for some r n related to the most relevant scaling operators. 



V. MULTIFRACTALITY AND ON-SITE PERTURBATION 

A. Density of states near the band center 

An important application of the analysis done in the previous sections is in the following. 
By expanding in E the action (j2J), composite operators like E n Tr(Q n ) appear in the model. 
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These operators can be written as linear combinations of scaling operators and among them, 
for each n, the most relevant one has dimension 

(52) 



z n = 2 + g\^- + 2jn - gn. 

Considering only the operator with this dimension as the most representative for E n , in the 
limit s — > oo, using Eq. (|50|). we can write 

"In s 



In 



~E n (s)~ 




E n 


-J 

Jo 



z„ dins'. 



Let us now define the following distribution function 

1 



□ 



yins 



/(In y/ Ins) 



(53) 



(54) 



where the function f(a) is defined by 

f(a)=2-(a + g) 2 /(2g(T + A)), (55) 

linked to z n by a Legendre transform (2^. Indeed finding a as solution of the following 
equation 

/'(a)+n = 0, (56) 

we obtain 

a(n) = ng(T + 4)-g. (57) 
This implies that z n and f(a) are related by the following Legendre transform 

z n = na(n) + f(a(n)). (58) 

Now we define through the distribution function (p>fj) the following mean value 

(y n ) ny ,s) = J dyv(y,s)y n . (59) 

Changing the integral variable by a = ln^/lns, Eq. (|59jl becomes 

(y n )v(y,s) ~ J das na+f ^ ~ s z ", (60) 

where we have evaluated the integral by the saddle-point method, being (jBTf) the saddle 
point. Using this result we can write E n (s) in terms of the average value of y n 

E n (s) 



E n 





" pins 






~ exp 


/ z n d In s' 


~ exp 


[/ 




yo 







" lns / dins' 



(61) 
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In some conditions this quantity is dominated by the tails of the distribution, namely it is 
determined by rare events and, therefore, can not represent the true energy scaling. From 
Eq. (jnOj) we notice, indeed, that for large n the tails of the distribution affect strongly the 
mean value of y n . The so called typical mean values, on the contrary, give insight on the 
bulk of the distribution function. For this reason we will consider the following quantity 



J typ\°) 

— — — exp 

^typ 



In s 



# p d\ns 



o 



exp 



In s 



o 



typ 



dins' 
In s' 



(62) 



where now z^ p is the typical dimension to be determined and (y n ) is replaced by y^ yp = 
exp (In y n ) with 

(\ny n ) H y )S) = J dyV(y,s)lny n ~ n \ns J daas^ a \ (63) 

where the integral can not be extended to — oo otherwise we would obtain irrelevant opera- 
tors. However from Eq. (J63)) we can see immediately that the regions of a where f(a) < 

E n (s) 

give negligible contributions to the typical value — ^ — . Let us consider the solution of 
/(«) = o 

a = 2 v^(f+4) - g = n c g{V + 4) - g, (64) 

discarding the other solution that is irrelevant in all regimes of disorder. In Eq. (J64|) we 
have introduced the factor 

n c = — . (65) 

in order to write at in analogy with Eq. ([57)1 . A more convenient way to remove from 
consideration all the rare events with large \a\ but with very small weig ht srt a ) with / (a) < 
is to write y™ yp in the same form of Eq. (|6()j) but with a restriction in the integration range 

/ das na+f{a) (66) 

Jf(a)>0 

since in the region of f(a) > we can expand s na ~ 1 + nalns. The same definition of 
typical values expressed by Eq. (jHEJ) has been already used in Ref. |l^ for the typical inverse 
participation ratios. 

We find that for n > n c the saddle-point (|57|) is outside the integration domain since 
f(a) < 0. The main contribution in the integral (jB^jl is then due to the boundary a, that 
do not depend on n, implying 

z y = na = n(2^g(T + A) - g). (67) 
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For n < n c instead the integral is determined again by the saddle-point (jBTfl since a is inside 
the integration domain being f(a) > 0. In this case we obtain z l ^ v = z n . 
Summarizing we have the following typical dimension 



ytyv 



2 + g + 2) n 2 — g n, for n < 



2nJg(T + 4) - an, for n > . 2 



The so called dynamical exponent |12] defined by 



z(s) = zT/n 



(68) 



in the strong coupling regime, where T > 9 < , has therefore the following scaling behavior 
for all value of n 

z(s) ~ 2y^f~ AgoVhiS. (69) 

Calling A typ = ^/E^ yp (s) the upper energy cut-off and C = ^ a positive constant, we have 
finally 



A 



typ 



E 



cxp 



Hip 



In s 



exp 



C(lns) 



^(s') dins' 

Now we can easily calculate the density of states p from its scaling equation 

dp 



dins' 



[z(s') - 2} p. 



(70) 



(71) 



Integrating over the scaling factor up to s = s(A typ ), we obtain 

P(s) 



In 



Po 



C(lns) 



21ns. 



_Li n 



From Eq. (fTUj) we have 

s = exp 

obtaining for the density of states the following behavior in energy 

A* 



(72) 



(73) 



P{E t yp) 



Hyp 

P0T^ ex P 



E 



i UP 



c 



1 In ^ ^ tyv 



E 



typ 



(74) 



We find that the density of states shows a weaker divergence than that obtained by Gade 
and Wegner 0,0] who found the exponent | on the logarithm. The final result (EH), on the 
contrary, is in perfect agreement with the density of states predicted in Refs. 
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B. On-site disorder 



In the presence of small on-site disorder or same-sublattice regular hopping, terms like 
c nr TiQ 2n appear in the theory |8j, responsible for chiral symmetry breaking of the two sub- 
lattice models. Following the same steps described previously, we have, in this case, 

(75) 
(76) 



2 + Ag[ - + 2)n z -2gn, 



related to the function 

f(a) = 2-(a + 2g) 2 /(8g(T + A)), 
through the following value of a 

ct(n) = Ang(T + A) - 2g (77) 
by the Legendre transform in Eq. ([58)1 . Defining again the typical value for c n , we get 

z tup = n(Ay/g{r + A)-2g), (78) 
meaning that the dynamic exponent is 

z(s) ~ A^/gf ~ 8g V\ns (79) 

and 



A 



typ 



Ctyp 



cxp 
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3#o 



(Ins 



exp 



2C(lns 



(80) 



The coefficient 2C in Eq. (jSUj) is twice the value which appears in Eq. (J70J). For this reason, 
considering the two on-site perturbations due either to a finite potential energy E or to an 
on-site disorder with strength c, if E ~ c, the crossover from the chiral symmetry to the 
standard one occurs first in the presence of the latter source of symmetry breaking. 



VI. CONCLUSIONS 



In this paper we have computed the anomalous scaling dimensions of an infinite family of 
operators in a non-linear a-model induced under RG by an on-site perturbation. We have 
applied this analysis to calculate the density of states near the chemical potential. 

The new result of the present work is to prove that the same expression for the density of 
states, already obtained through other approache s [rI 1^ which take advantage of the RG 
method proposed by Carpentier and Le Doussal [2lj, can be derived also within the more 
conventional non-linear a-model approach based on the replica method. 
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APPENDIX A 



In this appendix we present a more general version of Eq. 



Defining the operators 



0$» = A mQ ^m+iQ 4»-iQ A n Q, for m < n, (Al) 

0<£„ = A m Q A m+1 Q A ni - X Q A ni Q A X Q A n ^Q A n Q, for m > n, (A2) 

where A n are some symmetric or antisymmetric matrices in replica and frequency spaces, 
we obtain the following average value over fast modes 



i=i 



1+1 2^n Si +n + n 2 



—L 



i=l 



M 



JjTr'0? 



i=i 



i>j kj^i,j h=l lj=l 

N / N rii-i m—l 

2E II< E E « w 

i=l p=l Z=p 



(A3) 



where, as before, n = Ylf' n i an d the numbers ns i are defined by 

n^ — 1 m rii 



»* = e e n ^ 

m=l fc=l j=fe+l 



(A4) 



in which S^A^) is a sign, +1 if Aj is a symmetric operator and —1 if Aj antisymmetric. If 
all Aj are symmetric and equal, namely Aj = A = A 1 Vj, as in the particular case described 
by Eq. (J2EJ), where A is the identity in all the spaces, one has simply n s . = \riiirii — 1). 



APPENDIX B 

Here we consider, as an example, the case with n = 4. The polynomial in Eq. (|25|) then 
is the following 



Pa 



A^Tr'Q 4 + X^Tr'Q^Tr'Q + A^ 2) (Tr'Q 2 ) 2 + A^ ^Tr'Q^Tr'Q) 2 + (Tr'Q) 4 
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and we have from Eqs. 



the following /3-functions 



A (4) 



d\ W + q 

(4) y 



(16 + 8r) A (1) - 12A (2) - 16A (2) 

V 1 W (3,1) (2,2) 



/3, (2) = d\ (2) +g (10 + 8r) A (2) - 8A (1) - 16A <3) 
Pax, = d\ m + g \(8 + 8T) A W - 4A W - 4A (3) 

X {22) < 2 ' 2 ) L (2 ' 2) (4) (2 ' M) 



= uA + q 

A (2,1,1) 3 

(2,1.1) 



(6 + sr) A* 



/3,(4) = c?A (4) + o 

*A, 1,1,1,1) 3 

(1,1,1,1) k ' 



(2,1,1) 
(4) 



6A (2) -4A (2) -24A (4) 

(3,1) (2,2) (1,1,1,1) 



(4 + 8r)Ar -2A (3) 

V ' (1,1.1,1) (2,1 



Solving the four independent equations (J3U|) one can obtain 

the following real solutions (^ 4 > ^[3 1) ' ^ 22): ^21 1) ' ^[11 1 1)) = 
{ Cl (-48, 32, 12, -12, 1), c 2 (8, 4, -2, -5, 1), c 3 (2, -8, 7, -2, 1), c 4 (-4, -2, -2, 1, 1), c 5 (6, 8, 3, 6, 1)} 
with Cj arbitrary constants. Inserting these solutions in Eq. (}3*2"j) we obtain five scaling 
operators that behave in the following way 



(Ox) = (1 + (28 + 8T)L) 0i, 

(0 2 ) = (l + (14 + 8r)L) o 2 , 

(0 3 ) = (l + (8 + 8T)L) O s , 

(0 4 ) = (1 + (2 + 8T)L) C» 4 , 
(G 5 ) = (l + (-8 + 8r)L) e> 5 . 

For all of them Eq. (J3~3~j) is verified: Oi is related to the partition {4}, 2 to {3, 1}, O3 to 
{2, 2}, 4 to {2, 1, 1} and finally 5 is related to the partition {1, 1, 1, 1}. 

Alternatively from the /3-functions written above we can construct the matrix M which 
appears in Eq. (|3*T|) 



M 



( 


16 


-12 


-16 





o\ 




-8 


10 





-16 







-4 





8 


-4 










-6 


-4 


6 


-24 


V 











-2 





8n, 
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where I is the 5x5 identity matrix. M is diagonalized by the invertible matrix 



/ 


A O 

—48 


8 


2 


4 


6 1 




32 


4 


-8 


-2 


8 




12 


-2 


7 


-2 


3 




-12 


-5 


-2 


1 


6 




1 


1 


1 


1 





and its diagonal form is 

T~ X MT = diag{28, 14, 8, 2, -8} + 8H. 
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